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The classical formalism of the Moment Problem has been combined with a cumulant 
approach and applied to the extensive many-body problem. This has yielded many new 
exact results for many-body systems in the thermodynamic limit - for the ground state 
energy, for excited state gaps, and for arbitrary ground state averages. The method 
applies to any extensive Hamiltonian system, for any phase or symmetry arising in the 
model, whether on a lattice or in the continuum, and for any dimensionality. The 
theorems are of a nonperturbativc nature with respect to any couplings occuring in the 
model. 



1 The Early Development of the Lanczos Algorithm 

The Lanczos algorithm or method has been of interest to physicists because it is an 
essentially non-perturbative approach to physical problems with strong coupling, 
such as occur in the extensive many-body systems of condensed matter physics. 
In this method the Hamiltonian is used to generate a sequence of orthonormal 
states {\ipn)}n=i,2.. and Lanczos coefficients {an}n=o,i.., {Pn}n=i.,2.., from a suitably 
chosen trial state {ipo) through the following recurrence 



IV'n+l) = ^T— K-H" - an)\lpn) " f3n\tpn-l)] 
Pn+1 

SO that the Hamiltonian in this new basis is tridiagonal 
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In the mathematical and statistical arena the Lanczos Process has been long un- 
derstood as one manifestation of a bodjL of intimately connected mathematkial sub- 
jects, namely -rtiie Moment ProblcntJa, Orthogonal Pobjnpmial SystemiJo, Pade 
Approximationflu, Analytic Continued Fraction Theorjo'EI and Krylov Subspace 
MethodiO. One such equivalence is that the Lanczos process applied up to tit itera- 
tions to generate the Lanczos coefficients {an}"=0' {/5n}"=i precisely equivalent 
to generating the first Iut + 1 moments {p^r^Uv^o"^ defined by ^„ = (i?") ((O) 
denotes the expectation value with respect to the trial state). 

The traditional use of the Lanczos algorithm has been in a purely numerical 
way, that is to say as a numerical technique for exact diagonalisation of very iarge 
matrices that arise in treating many-body problems in small finite systemfl, or 
in the treatment of the one-electron problem in disordered or aperiodic systems, 
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as in the Recursion Method of Haydoclfl. The potential of taking the Lanczos 
algorithm far beyond these limitations, into a more powerful, universal formalism 
has not been widely jappreciated, although some inkling of this was apparent in the 
suggestion of Mattifj concerning the exact mapping of the many-body problem 
onto a one-dimensional nearest-neighbour model. This idea was £i!a|ilD|md in some 
applications to the Kondo and Wolff models by Mancini and MattiO'E5E^. We wish 
to emphasis to the reader that our approach here is quite different from that used 
in the exact diagonalisation studies of finite systems in two respects - we do not 
construct a full basis for a finite system but manipulate basis vectors and coefficients 
of an arbitrarily large system analytically and symbolically, and we perform every 
iteration exactly and therefore need not concern ourselves with round-off or loss of 
orthogonality issues. 

While little use or development of the mathematical constructs were employed 
in exact diagonalisation methods, some of the ideas were used in other £amW|lispis. 
Formalisms weie, desieloped under the name of the Recursion MethooB'liallalla or 
related methodji3't§llj for one-electron problems, but in the last analysis every cal- 
culation was a numerical evaluation, i.e. explicit construction of the orthogonal 
polynomials via 3-term recurrences and then the continued fraction representation 
of the density of states. Some questions were raised concerning the generalisation 
to genuine many-body problem£3 but this was not realised at the time. Other 
formalisms, treating many-body systems and stochastic processes in the thermody- 
namic limit, ..which arose from this mathematical legacy were the Memory Function 
formahsmEJc3€3'oE3, the Recursion MethocO (not iflpbe-caDfLwed with the previ- 
ous use of the same term) and the Projection MethodEZln3'E3E3E2l. However all these 
methods were applied only formally, that is to say the consequences of introducing 
these tools into the many-body problem was not systematically followed through 
or explored - the recursion process would be carried out up to a finite number of 
steps and truncated in an ad-hoc manner. This can be done analytically by hand 
for the first few steps, but usually higher steps are calculated on a computer by 
constructing an equivalent graphical description of the problem and making the 
combinatorial evaluations that arise. The formalisms are precise and exact in this 
regard but being truncated in this manner they have not converged ut —>■ oo. This 
is a serious issue because while nr may be numerically large, say 20 or 30, one also 
wants to follow this Lanczos convergence with the thermodynamic limit iV — + oo, 
but the problem is that the value of Uopt, to assure convergence to a given accuracy, 
will scale with N at best, and may sometimes scale with a higher dependence. 



2 The New Developments 

However it is possible to transcend these limitations in the process of constructing 
the mathematical formalism properly embedded in its physical context. The first 
key ingredient is to find a way of incorporating the system size scaling for the 
extensive system into the existing formalism from the outset. The solution to this 
is obvious-^i describe everything in terms of cumulants, connected moments or semi- 
invariant£3 {j^njJJ^i {i^n = {H"-)c) instead of moments. The defining relationship 
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and there exists a unique transformation between the set of first ut cumulants and 
the set of first jit moments. Cumulants scale with the system size in the following 
way 

l^n = CnN , l^n = C„iV + m„ , (4) 

in the ground state sector and other sectors respectively, ignoring boundary con- 
ditions. The coefficients c„, to„ are independent of N and functions of coupling 
constants and other parameters in the trial state. With this scaling form all infor- 
mation regarding finite-size scaling is lost, but it is the simplest approach. Uncon- 
nected moments encapsulate the information about a system in a very redundant 
way and which leads to problems of ill-conditioning. 

Once the above step is taken then many results become quickly apparent. The 
first result arises from the substitution of the cumulant Eq. into the explicit 
expression for the moments, and then into the Lanczos coefficients, and expanding 

larkably simple 
] - as a function 



the resulting forms in a large N expansion - and the result is_a 
and perfectly universal expansion - the "Plaquette Expansion'^ 
of an arbitrary Lanczos iteration number n 
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It can be showrE^ that just retaining the first terms in each coefficient can be in- 
terpreted as a manifestation of the Central Limit Theorem and a description in 
terms of Gaussian fiuctuations, while retaining the first two terms in each describes 
dynamical processes governed by the Binomial distribution and all related ones. In 
this way many of the distributions of statistics arise naturally in such an expan- 
sion, and in fact simple interacting physical models can b£_found which are exactly 
represented by a finite number of terms in this expansiortj. 

It should be noted that in each term of the above the degree of the polynomial 
in n is the same as the inverse power of N so that the following limit n, N ^ co 
exists at fixed s = n/N. Although the above is just a Taylor series expansion in 
1/A'' we conjecture that the exact Lanczos coefficients exhibit the following general 
confiuence 

a^N) -^-^^ Na{s) , (3l{N) "'^-7°^ N'(3\s) . (7) 
In the next step if one defines the spectral envelope functions 
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then one can employ theorems on bounds to the extremal zeros oLiDrthogonal 
Polynomials to arrive at an exact theorem for the ground state energjEj 



eo = inf [e(s)] 



(9) 



and if this occurs at a finite point it is denoted sq. This result constitutes an 
exact diagonalisation of the many-body problem in the thermodynamic limit, as the 
formalism expresses results in terms of the tridiagonal matrix elements, or Lanczos 
coefficients. From this, it is a simple step to finding the ground state average for an 
arbitrary operator ($3 
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where the operator Lanczos coefficients are constructed from the operator cumulants 



-1 = ^{H''-^dH^)c S"a,S°P^ 



(11) 



fc=0 



in the same manner as the pure Lanczos coefficients. The excited state gap, between 
the ground state and an excited state in another sector, is just tke difference between 
two ground state energy densities, ei — eq = Ae/N and is thufia 



P{s) 



(12) 



where the gap Lanczos coefficients are constructed from the gap cumulants (con- 
structed using a trial state with the excited state quantum numbers) 

i^n = CnN + 5%, ^ S'^a, S^P^ . (13) 

For the excited state gap in the same sector, the following peeling theorem holdfS 

Ae = 2 hm iV [e(s) - e„(iV)],^ . (14) 

In its application to non-integrable models the above expansion, Eq.(^,^, is gen- 
erated from a finite set of low order cumulants and then truncated at some finite 
order and the above theorems applied without the need for any extrapolation. Some 



examples where t 
Heisenberg models 



las been successfully enmlajzied are the 1 and 2-dimensional 
I and lattice gauge modclMSH. There are also examples of this 
method used in an essentially exact manner namely for a 1-dimensional solvable 
spin model with a phase transition at T = (o, where the convergence properties of 
the method have been examined. 



3 Summary 

As we have seen the Analytic Lanczos Method is an important stage in the develop- 
ment of Lanczos methods in the extensive many-body problem. Amongst its virtues 
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are that is general purpose - it works for any Hamiltonian, lattice or continuum, 

quantum mechanical or classical, in all dimensions D, for any phase or symmetry 
of the model, it is non-perturbative in couplings, it works exactly in the thermo- 
dynamic limit N ^ oo and it applies to ground state or T > properties. It is 
accurate and systematic in that there is a development in successive orders so that 
some control of the errors can be made. It has a flexible implementation in that 
the treatment can be either analytic, semi-analytic or numerical depending on the 
degree of integrability of the model at hand, that one is free to choose the trial 
state, within very general limits relating to the symmetry of the target state, and 
that one can combine it with other methods, e.g. variational, mean-field, 
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